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^ ' Abstract 

(N ■ 

^ ' We propose bigradient phase referencing (BPR), a new radio-observation technique, and report its 

^>0 , performance using the Japanese very-long-bascline-interferometry network ( JVN) . In this method, a weak 

source is detected by phase-referencing using a primary cahbrator, in order to play a role as a secondary 
calibrator for phase-referencing to a weak target. We will be given the opportunity to select a calibrator 
from lots of milli-Jansky sources, one of which may be located at the position closer to the target. With 
1^5 ! such a smaller separation, high-quality phase-referencing can be achieved. Furthermore, a subsequent 

f^ . more-sophisticated calibration can relocate array's focus to a hypothetical point much closer to the target; 

a higher quality of phase referencing is available. Our demonstrative observations with strong radio sources 
^ , have proved the capabilities of BPR in terms of image dynamic ranges and astrometric reproducibility. 

I ' The image dynamic range on a target has been improved with a factor of about six compared to that of 

O I normal phase-referencing; the resultant position difference of target's emission between two epochs was 

4— » ' only 62±50 micro-arcsecond, even with less than 2300-km baselines at 8.4 GHz and fast-switching of a 

~^ \ target-calibrator pair of a 2.1-degree separation. 

• • , Key words: astrometry — atmospheric effects — phase-referencing — techniques: interferometric — 

J^ ' very-long-baseline interferometry 
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1. Introduction The first solution can be achieved by geodetic- like obser- 
vations (Brunthalcr, Rcid, & Falcke 2003; Mioduszewski 
A phase-referencing technique allows very-long-baseline & Kogan 2004; Reid & Brunthaler 2004) or by parallel 
interferometry (VLBI) to make relative astrometry with plate air modeling of the long-term phase drifts (Reid et 
an accuracy of less than 1 milli-arcsecond (mas) and to al. 1999; Brunthaler, Reid, & Falcke 2003). The data pro- 
detect very weak sources in mJy level. The quality of cessing of the geodetic-like observations is supported by 
phase referencing is limited by residual errors in differen- the task DELZN of the Astronomical Image Processing 
tial phases between a target and a calibrator (e.g., Bcasley System (AIPS; Greisen 2003), developed at the National 
& Conway 1995). The most important error compo- Radio Astronomy Observatory. The second solution can 
nent is the uncertainty of atmospheric models in corre- be achieved by using several strong calibrators over sev- 
lators. The differential excess-path length between the eral degrees; this calibration process is also supported by 
two sources at different elevations is significantly harm- the task ATMCA in the AIPS (Fomalont & Kogan 2005). 
ful, even with a zenith phase-delay error of only a few The third solution, using a calibrator very close to a 
centimeters and a separation of sources only one degree target, depends on a matter of blind chance. However, it 
(Beasley & Conway 1995; Reid ct al. 1999). Possible so- is promising if array's sensitivity is significantly improved, 
lutions are (1) estimating the unknown phase-delay at since the surface number density of radio sources dramat- 
zenith, (2) to determine the residual phase gradient in ically increases in the fainter sky (e.g., Fomalont et al. 
the sky, and (3) using a calibrator very close to a target. 1991). The present paper proposes a bigradient phase- 
It is advisable to apply all these solutions simultaneously, referencing (BPR) method, which gives us a chance to 
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utilize such weak radio sources as calibrators. Once a 
weak source very close to a target is detected by phase- 
referencing using a strong calibrator, and then plays a 
role as a calibrator for phase-referencing to detect a very 
weak target. Such a two-step approach is responsible 
for its name. Furthermore, assuming that the residual 
phase components are linear around the target (Fomalont 
& Kopeikin 2003), array's focus on the weak calibrator 
can relocate to a point much closer to the target. Hence, 
the BPR leads to be nearly free from the long-term phase 
drifts. 

In the present paper, we present the BPR method and 
its observational tests. In Section 2 wc describe the prin- 
ciple of the method and predict final phase errors. In 
Section 3 we report demonstrative VLBI observations. In 
Section 4 we discuss the capability of this method. Finally, 
we summarize the method and the test observations in 
Section 5. 

2. Method 

2.1. A calibrator arrangement and scheduling 

The BPR method schedules three sources: a target 'T', 
a primary calibrator 'CI', and a secondary calibrator 'C2'. 
CI is bright enough to detect in half of a coherence time, 
typically a few minutes or less, depending on observing 
frequency and weather condition (Ulvestad 1998). Wc as- 
sume here C2 provides correlated flux densities insufScient 
for fringe detection in half of the coherence time but suf- 
ficient for detection by phase-referencing of several tens 
of minutes: it is ~10 mJy at centimeter bands for typ- 
ical VLBI arrays. On the basis of statistical discussion 
on the surface-number density of radio sources in the sky 
(e.g., Fomalont et al. 1991), CI may be located relatively 
far from T, typically ~2 deg separation, while C2 can be 
found at the position significantly closer to T, possibly 
<G.5 deg separation. An example of source configuration 
is shown in Fig. 1. 
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Fig. 1. An example of source configuration on the sky 
plane, and scheduling. Only a primary calibrator CI is 
fringe-detectable. C2 is a secondary calibrator, a weak source. 
C2' is a hypothetical calibrator, on blank sky at the tangential 
point on C1-C2 line. 

Both T and C2 should be observed using phase- 
referencing because of their weakness. An example of 
schedule for the BPR is as follows: 



. . . - CI - C2 - CI - C2 - CI - C2 - CI - C2 
-CI -T-Cl -T-Cl-T-Cl -T 
-CI - C2 - CI - C2 - CI - C2 - CI - C2 
-CI -T-Cl -T-Cl -T-Cl -T-... (1) 

It consists of sets of fast-switching for two pairs: C2-C1 
and T-Cl. The schedule intends to remove rapid phase 
fluctuations by fast-switching, and to change the reference 
point CI into C2 by pair-swapping in order to reduce the 
separation angle from T. The scans for the C2-C1 pair 
must be allocated every less than 1 hour, so that the long- 
term phase drifts can be tracked. The number of iterations 
of fast-switching for the C2-C1 pair should be set to secure 
the detection of C2 with a signal-to-noise ratio of more 
than 5, for successful self-calibration in phase domain. 

2.2. Observation equations 

With such an observing schedule, an observer will ob- 
tain raw visibility data for the target and two calibrators. 
After determination of amplitude-gain, delay and delay- 
rate solutions using CI, all we have to consider is phase 
terms. Observed, raw visibility phases, (t>oBS, involve var- 
ious phase terms. 

(2) 
(3) 
(4) 
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(fp- are phase terms for CI, C2 
'stru is a phase term originating 



where 

and T, respectively, 
in source structure; (j)pos is a positional-phase delay rel- 
ative to a phase-tracking center; (/)inst is an instrumental- 
phase delay; (j)gco is a geometric-phase delay error; (patmo 
is tropospheric/ionospheric-phase delay error; (/irapid is a 
rapidly variable phase due to water vapor flowing at low 
altitude. We ignore thermal phase noise and calibration 
errors in the present paper. We deal with the phase terms 
of source position and structure separately, in order to 
show an astrometric term and its error expressly at the 
last equation. Time scales of the change of these phase- 
error components are usually more than several tens of 
minutes, except for (/)rapid in which the time scale is typi- 
cally a few minutes or less. 

2.3. Calibrations of bigradient phase referencing 

Data of CI have high signal-to- noise ratios enough to 
do self-calibration. Hence, a sufficiently feasible source 
structure model can be obtained: 
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(5) 



(6) 



An antenna-based solution table, provided from the self- 
calibration, includes all the terms other than the source 
structure. This is a 'solution table 1,' ^sNi- We express 
a determined term in its capital, in the present paper. 

With this solution table, we can do phase correction to 
both C2 and T, so that their phases are free from rapid 
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Table 1. Array configurations for our observations. 



Epoch Date 



Telescopes 



1st 2005 Sep 25 08:00-12:30UT 

2nd 2005 Oct 23 06:00-12:00UT 



VERA X 4 (Mizsawa20ni, Ogasawara20m, Iriki20m, Ishigaki20m) 
VERA X 4, Kashima34ni, Tsukuba32m. Usuda64m, Yamaguchi32m 



phase fluctuations due to water vapor. The calibration on 
C2 is as follows, 
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where A(/)p-ci ^ ^p 

/C2-C1 



+ A0,^,;Tj"\ (7) 

>f^, differential phase terms. 
Presumably, A0j:^lT^^ « 0, because its random fluctuation 
is not responsible for a systematic residual, so should be 
averaged out. Also Ac/)?^^^^ w 0, because an identical re- 
ceiving system is used. The calibration on T can be done 
in the same manner; this is a normal phase- referencing, 
which has been commonly used (Beasley & Conway 1995). 
The other differential phase terms may not be zero. 
They are responsible for long-term phase drifts, which 
may cause some apparent position shift and the degra- 
dation of image dynamic range. As we mentioned above, 
the time scales of change of these errors are usually more 
than several tens of minutes; we can integrate the data co- 
herently in a period of less than their time scales. When 
signal-to-noise ratios sufficient to do self-calibration are 
available on C2 by this integration, its source structure 
model can be determined. Then, equation (7) becomes 

-*SN2, (8) 
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Thus, we obtain a 'solution table 2,' <i>sN2- We recom- 
mend that the first self-calibration on C2 is done using 
a tentative point source model at the phase-tracking cen- 
ter, in order to include (pp^^ into the solution table 2 for 
astromctry. 

Now we apply both the solution table 1 and 2 to the 
target data; 



'OBS 


- ($SN1 + $SN2) 




= (fltTU + (t^pos + '/'hist + '/'Jeo + '/'Itmo + '/'rapid 




- ['/'pos + '/'Pnst + '/'gco + </>Smo + '/'rapid 




+A'^p°o\-^^ + A0^c-c^-fA0S^„-c^] 




= '/'In + A^^-c^ + A0j-c^ + A0j,-Si 




+[A<^T;C^ + A0l,-?]. (10) 



This equation is very similar to equation (7), but this in- 
cludes the T-C2 pair. Again A(j)l~^^^ « and A(j>J^-^^^ « 
0. The process in this equation has replaced CI with C2 



as a reference point, where array's focus is strictly opti- 
mized. In other words, we can make a substantially phase- 
referencing for the T-C2 pair without fast-switching be- 
tween them. A fringe-undetectable source becomes use- 
ful as a reference calibrator. The separation angle of 
the T~C2 pair is smaller than that of the T-Cl pair, so 
A0J-C2 < A0T-C1 and A^Jj-'^i^ ^ A^Jj-'^i^ This means 
that the differential excess path becomes smaller, and we 
would obtain a better image quality in comparison with 
that of the normal phase-referenced image. This is the bi- 
gradient phase-referencing. The astrometric measurement 
of T can be done relative to C2, because of A^p"'-'^, with 
the uncertainty of Ac/iJ^"^^ -I- AcjyJ^^^ , as shown equa- 
tion (10). 

In most cases, observers have to be careful about (j)^^^ 
in the stage of self-calibration on C2 [equation (8) and 
(9)]. The mas-scale position of C2 may be unknown 
when scheduling, because such a weak source is not cat- 
aloged into the International Celestial Reference Frame 
(ICRF; Ma et al. 1998). In the phase-referenced image of 
C2 [equation (7)], an emission peak may initially appear 
shifted by several tens of mas, even if the phase-tracking 
center was set based on a position measured with the Very 
Large Array (VLA) A-array configuration, etc. This leads 
a too large $sn2 to make phase-connections beyond 2tt 
ambiguity. Also, the requirement of position accuracy of 
C2 is the same as that of CI, because C2 has become a 
reference calibrator as a substitute for CI in the BPR. 
The large position error of a reference calibrator will de- 
grade the dynamic range of target image (section 17.3.5 
of Beasley & Conway 1995). We recommend correcting 
the phase-tracking center to an accurate position by post 
processing with the tasks CLCOR or UVFIX in the AIPS, 
before the self-calibration on C2. 

2.4- A more- sophisticated calibration 

The solution table 2, <E'sn2, has tracked the long-term 
phase drift between the C2-C1 pair, shown as equa- 
tion (9). However, we hope ultimately to know that be- 
tween the T-Cl pair. The method of Fomalont & Kogan 
(2005) achieves it using more than one calibrator around 
a target, in the assumption that long-term phase varia- 
tions arc linear over the region of the sources (Fomalont 
& Kopcikin 2003). We also apply the same assumption 
to an additional calibration, which is more-sophisticated 
than the BPR. Array's focus has moved to C2 from CI by 
the BPR (section 2.3). This means the solution table 2 
has the ability to shift the focus. We can make a new focus 
at any point on the line of C2-C1 as far as the linearity is 
valid. In the more-sophisticated calibration, we establish 
the hypothetical calibrator C2', on blank sky at the tan- 
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gential point on the C1-C2 line from T, by scaling $sn2 
with a factor r, vfhere r = [O'^^^'^'^' / O'^^^^'^) , 0Ci-C2' ^^-^^ 
^Ci-C2 gj.^ separation angles of C1-C2' and C1-C2, re- 
spectively (Fig. 1). If both the solution table 1 and this 
modified solution table 2 are applied to the target data, 
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Again A(/)jy^7j^^ «0 and A^^'^;^*^^ wO. Because of the linear- 
ity, the terms r x A^g^^^^^ and r x A(l)'^^~^^ virtually cor- 
respond to differential geometric and atmospheric phase- 
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. Finally, equation (11) becomes 

r X $SN2) 

r X A<„^-Ci 

geo • A0T-C2' ], (12) 

This is substantially a phase-referencing for the T-C2' 
pair, although there is no fast-switching between them. 
The modified phase solutions of self-calibration have made 
array's focus at the sky toward C2', much closer to the 
target. 

The astrometric term in equation (12) is A(/)p~^^. The 
term r x A(j)p^~'-^^ consists of only positional terms, so 
it makes only a slight position shift without degrading 
an image dynamic range. Even with ICRF calibrators, 
their absolute-position uncertainties ('^ 0.3 mas, Ma et 
al. 1998; Fey et al. 2004) may be responsible for the term. 
Since the differences between the real positions and phase- 
tracking center are unpredictable, r x A^p^"*-^^ is also un- 
predictable. However, this will practically not contribute 
to astrometric errors in observations of relative astrom- 
etry, because r x A(j)p^~'^^ is duphcated every epoch if 
the relative positions between the emission peaks of cal- 
ibrators and the phase-tracking centers are not changed 
any epochs. Thus, we can make relative astrometry for 
measuring proper motion in series of monitoring observa- 
tions. Please note that observer must use identical phase- 
tracking centers every epoch and the first self-calibrations 
using the structure models at the phase-tracking centers. 

The terms of [A<l)^^'~'^ + Acj)^^^-^^^ ] are expected to be 
zero, if the position of C2' is coincident with that of T, 
i.e., all the three sources align exactly on a straight line. 
If the sources are misaligned, these terms become a phase 
error, depending on the separation angle between C2' and 
T. Because the BPR makes a lot of weaker calibrators 
usable, observers will have a desirable allocation of sources 
in the sky much easier. Possible solutions to reduce this 



remaining error of the misaligned case are discussed in 
Section 4.3. 

2.5. Calibration errors 

In addition to the errors theoretically expected in the 
method, there are various calibration errors. These kinds 
of errors are usually not so large compared to the long- 
term phase drifts, but still sufhcient to affect the quality 
of a phase-referenced image. Observers must be careful to 
address them as much as possible, when pursuing accurate 
astrometry. 

(1) The purity of the solution tables is related to the 
quality of self-calibration in equations (5) and (8). Very 
complex source structures or low signal-to-noise ratios 
make it difhcult to separate $stru from the other phase 
terms. This calibration error would be significant, because 
the BPR involves self-calibration twice, and because one of 
two calibrators is always a weak source. An antenna-based 
phase solution with a signal-to-noise ratio of 3 corresponds 
to an accuracy of ~ 21 deg theoretically (Thompson et al. 
2001), which would be as large as systematic phase com- 
ponents to be corrected. Hence, we recommend the use 
of C2 with an signal-to- noise ratio of more than 5. Note 
that phase signals are usually weaker on longer baselines. 

(2) The rapid phase fluctuations with the same as or 
less than the period of fast-switching cannot be removed, 
although we assumed Acj}'^^^^^'^ w and AcjjJ^^^^^ « 0. The 
remaining fluctuations tend to be random and averaged 
out, but are responsible for coherence loss. This situation 
also arises in normal phase-referencing of fast-switching. 
A dual-beam observation, as the VERA, is presumably 
the most effective ways against such a random phase fluc- 
tuation. 



(3) The longer term changes of 



■^pos; V^atmoi 



and 



have time scales of more than several tens of minutes. 
The pair-swapping tracks them, but cannot monitor their 
fluctuations with less time scales due to absence of fast- 
switching between T and C2. The fluctuations affect not 
only phase measurements in C2 scans but also calibra- 
tion parameters that will be made for T by interpola- 
tion between the C2 scans. We recommend that the pair- 
swapping period is about 20 minutes, which provides some 
redundancy in the time scale of long-term phase drifts, 
when accurate astrometry is required. The fluctuations 
will be averaged out, although some cohelence loss is not 
avoidable. 

(4) The structure change of calibrators results in a posi- 
tion shift of the target; astrometric errors of 100 microarc- 
seconds over observations longer than several months are 
possible. Although this is not a kind of calibration errors, 
it is not avoidable for observers. The intrinsic position 
change on a calibrator will bring the same position shift 
on the target without degrading its image dynamic range. 
It is necessary to use calibrators with stable structures, 
preferably point sources; do not use super-luminul quasars 
even if high signal-to-noise ratios are expected. 
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3. Test observations 

Wc present demonstrative VLBI observations by the 
BPR method in this section. Although this method aims 
to detect a weak target using a weak caHbrator and a 
strong caUbrator, for easy inspection of visibiHty-phase 
correction, here we selected a target and two calibrators 
from strong ICRF radio sources, as follows: 

• 3C 345 as a target 'T' 

• NRAO 512 as a secondary calibrator 'C2' 

• DA 426 as a primary calibrator 'CI.' 

The relative locations of these sources are similar to Fig. 1, 
but align more straightly on east-west direction on the 
sky plane. CI and T, a primary pair, are separated by 
2.09 deg; CI and C2, a secondary pair, are separated by 
2.57 deg; C2 and T are separated by 0.48 deg. T and C2' 
are separated by only 2.1 arcmin. 

3.1. Observations 

VLBI imaging observations with the BPR method were 
carried out at two epochs with the Japanese VLBI net- 
work (JVN), a newly-established VLBI network with ~ 50- 
2300 km baselines across the Japanese islands (Fujisawa 
et al. in prep). This array consists of ten antennas, in- 
cluding four radio telescopes of the VLBI Exploration of 
Radio Astrometry project (VERA; Kobayashi et al. 2003). 
The telescope participants for the observations are listed 
in Table 1. Right-circular polarization was received at X- 
band. Two frequency bands, 8400-8416 MHz (IFl) and 
8432-8448 MHz (IF2), were selected. The VSOP-terminal 
system was used as a digital back-end; digitized data in 
2-bit quantization were recorded onto magnetic tapes at 
a data rate of 128 Mbps. Correlation processing was 
performed with the Mitaka FX correlator (Shibata et al. 
1998) at the National Astronomical Observatory Japan. 

The switching-cycle period of phase-referencing was 5 
minutes. A set of switching cycles for the primary pair (T- 
CI) spent about 20 minutes; a set for the secondary 
pair (C2-C1) also spent about 20 minutes. The two sets 
were scheduled alternately, as mentioned in Section 2.1. 
Thus, a pair-swapping cycle period was about 40 minutes. 

3.2. Data reduction 

Data reduction procedures were performed for three 
data sets: Ist-epoch data (VERA antennas), 2nd-epoch 
data from only the VERA antennas, and 2nd-epoch data 
from all the antennas (Table 1). The correlated data 
were reduced using the AIPS. After initial data inspec- 
tion and flagging, fringe-fitting was performed to 3C 345 
and NRAO 512. Note that in this stage we did not in- 
tend to use the two sources as T or C2 but for bandpass 
calibration and amplitude-gain calibration, respectively. 
Since NRAO 512 is the most ideal point source up to 
our longest baseline (^ 50 MA), self-calibration in terms 
of amplitude provided us the solution of amplitude-gain 
variation relative to each antenna with the accuracy to 
about 1%. We have done absolute- amplitude scaling us- 
ing NRAO 512 with 0.89±0.09 Jy, a single-dish measure- 



ment with the Yamaguchi 32-m radio telescope (Fujisawa 
et al. 2002) at 8.38 GHz on 2005 Nobember 12. After 
these processes, we deleted the solutions of fringe-fitting 
for 3C 345 and NRAO 512. Data obtained at elevations 
of less than ~ 25 deg were flagged out. We utilized no a- 
priori gain parameter in both amplitude and phase except 
for the absolute flux scaling. 

Next, we performed fringe-fitting for CI (DA 426), then 
established its source structure model in the DIFMAP 
software (Shepherd 1997) using deconvolution and self- 
calibration algorithms iteratively. With the CALIB task 
in the AIPS, the phase solution of self-calibration was ob- 
tained using the source structure model: we obtained a 
'solution table 1' (section 2.3). We applied this table to 
the data of both T and C2. At this stage, we can obtain 
normal phase-referenced images. 

The phase-referenced visibilities of C2 (NRAO 512) 
were coherent in at least several tens of minutes (described 
in Section 3.3.2), as we have assumed. C2 was supposed 
to have sufficient signal-to-noise ratios to perform self- 
calibration in phase when integrated for several tens of 
minutes, even if with down to '^10 mJy. We have suc- 
cessfully obtained a 'solution table 2' by self-calibration 
for NRAO 512 in a solution interval of 40 minutes. The 
target image, which were synthesized from visibilities cor- 
rected with both the calibration table 1 and 2. are shown 
in Fig. 2 (b), (f) and (j); these are the images provided by 
the BPR method. 

Since the three sources quasi-perfectly align on a 
straight line on the sky. This is an ideal situation for 
the more-sophisticated calibration (Section 2.4). The 
phase solutions in the solution table 2 was scaled down 
to r (= [fl^^-^^'/^^^"^^]) times, and then the modified 
solution table was applied to the target data together with 
the solution table 1. The resultant images are shown in 
Fig. 2 (c), (g), and (k). For comparison, images corrected 
by self-calibration both in amplitude and phase are also 
shown in Fig. 2 (d), (h), and (1). 

The number of usable visibilities in the data set of 2nd- 
epoch from only the VERA antennas were ~ 30% less 
than the Ist-epoch ones, because the significant fraction 
of Cl's scans were flagged out in fringe-fitting. With such 
a small number of visibilities and a poor UV-coverage, 
the resultant images could not bear comparison with the 
other data sets. 

3.3. Results 

3.3.1. Image quality 

Images by the BPR appeared in Fig. 2 are dramatically 
improved compared to those by normal phase-referencing: 
image dynamic ranges have increased about 3.1-3.8 times 
(Table 2). The jet structure of 3C 345 is seen in the im- 
proved images. This is the fruit of using a closer calibra- 
tor adopted in the BPR. The normal phase-referencing 
was from the T-Cl pair, separated by 2.09 deg, while 
the BPR was practically from the T-C2 pair, separated 
by only 0.48 deg. a 4.3-times smaller separation. Better 
images have been obtained by the more-sophisticated cal- 
ibration. This is the fruit of optimization of the solution 
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Fig. 2. Contour maps of target 3C 345. (upper panels) Images at the first epoch with four telescopes of the VERA, (a) normal 
phase-referenced image, (b) Bigradient phase-referenced image, (c) Image made by more-sophisticated calibration. Contours of 
(a)-(c) start from SiTc, where CTc is RMS of image noise on the image (c) (see Table 2). (d) Self-calibrated image. Contours of 
(d) start from SfTdi where a^ is RMS of image noise on the image (d). Synthesized beams, shown at the lower- left of each image, 
are of 3.75 x 2.02 mas at a position angle of —69.3 deg. (middle panels) Images at the second epoch with only four telescopes 
of the VERA, (e) normal phase-referenced image, (f) Bigradient phase-referenced image, (g) Image made by more-sophisticated 
calibration. Contours of (e)-(g) start from 3(Tg. (h) Self-calibrated image. Contours of (h) start from 3(tjj. Synthesized beams are 
of 2.00 X 3.52 mas at a position angle of —67.3 deg. (lower panels) Images at the second epoch with all telescopes, (i) normal 
phase-referenced image, (j) Bigradient phase-referenced image, (k) Image made by more-sophisticated calibration. Contours of 
(i)-(k) start from 3(71;. (1) Self-calibrated image. Contours of (1) start from 'ia\. Synthesized beams are of 2.72 X 5.54 mas at a 
position angle of —70.6 deg. Contour levels for phase-referenced images are Scr x (—1, 1, 2, 4, 6, 8, 10, 12, 14, 16); for self-calibrated 
images are 3(7 x (—1, 1, 2, 4, 8, 16 ,32 ,64 ,128 ,256). Negative contours are described in dashed-curves. 



table 2 for the sky point of C2', separated from T by 
only 2.1 arcniin, a 61-tinies smaller separation than that 
of T-Cl. In the VERA image at the 2nd epoch, the more- 
sophisticated calibration appears not to work well in terms 
of a image-noise level. Because of small number of effec- 
tive visibility of this data set, measurements for this data 
set might be some inadequate. 
3.3.2. Long-term phase drifts 

In the BPR method, phase-referenced visibilities of C2 
are supposed to be coherent in at least several tens of 
minutes. Additionally, the more-sophisticated calibration 
expects that the phase of phase-referenced C2 is a linear 
extrapolation of the phase-drift between T and CI. We 
show an example of phase-referenced visibility phases both 
of 3C 345 (T) and NRAO 512 (C2) in Fig. 3. A long-term 
phase-drift is seen. The phase-drift of 3C 345 was undesir- 



able for normal phase-referencing. From DA 426's (Cl's) 
point of view, both the two sources are at about the same 
direction and separation. The phases of NRAO 512 are 
supposed to be similar to the drift of 3C 345, and this 
expectation met the obserevation. Consequently, we can 
remove a large part of the phase drift on 3C 345 by shifting 
array's focus to NRAO 512. 

The phase drift of NRAO 512 became larger than that 
of 3C 345 especially later than 09h, where elevations of 
these sources became lower. This deviation occurred be- 
cause that the separation angle of NRAO 512-DA 426 was 
larger than that of 3C 345-DA 426. Consequently, the 
data of 3C 345 will be over-corrected by the solution ta- 
ble 2 in the BPR. The more-sophisticated calibration aims 
to correct such a difference by optimizing the solution ta- 
ble 2 toward the hypothetical calibrator C2'. The RMS 



No. 



Bigradient phase referencing 

Table 2. Phase-referencing qualities in images of 3C 345. 



Epoch Ant. 



(1) 



(2) 



Method 



(3) 



(Jy beam^^) 
(4) 



Cimg 

(mJy beam" 

(5) 



/P/^in 

(6) 



0"^ 

(deg) 
(7) 



(mas) 
(8) 



A<5 

(mas) 

(9) 



1st 



2nd 



VERA 



VERA 



all 



PR 


3.49 


291 


12 


65 


-1.103±0.118 


+0.925 ±0.059 


BPR 


4.03 


107 


37 


26 


-0.956 ±0.054 


+0.238 ±0.025 


BPR+a 


4.06 


97.4 


41 


30 


-0.999 ±0.050 


+0.374 + 0.023 


self-cal. 


4.32 


9.23 


468 


21 






PR 


3.62 


290 


12 


69 


-1.308 ±0.159 


+1.198 + 0.074 


BPR 


3.90 


88.7 


44 


21 


-0.725 ±0.044 


+0.077 + 0.022 


BPR+a 


3.93 


100 


39 


28 


-0.871 ±0.050 


+0.307 + 0.025 


self-cal. 


4.10 


15.8 


259 


19 






PR 


3.04 


400 


8 


79 


-0.926 ±0.313 


+1.243 + 0.175 


BPR 


4.34 


144 


30 


33 


-1.001 ±0.083 


+0.126 + 0.043 


BPR+a 


4.56 


93.2 


49 


30 


-0.967 ±0.055 


+0.320 + 0.028 


self-cal. 


4.91 


7.26 


676 


18 







Comments — (1) observing epoch, (2) antennas used for mapping (3) calibration method. 'PR' represents normal phase referencing. 
'BPR' represents bigradient phase referencing. 'BPR-fo' represents more-sophisticated calibration, 'self-cal.' represents self-calibration. 
(4) peak intensity, (5) RMS of image noise, (6) image dynamic range, (7) RMS of visibility-phase scatter, (8) and (9) position shift from 
phase-tracking center of 3C 345. This was measured by two-dimensional Gaussian fitting using task JMFIT of AIPS. 



of calibrated phases are listed in column (7) of Tabic 2. 
Steady improvements were seen in the images from the 
2nd epoch data from all the antennas. However, slightly 
increases of RMS were seen in VERA data. It may be a 
coincidental effect of the intrinsic structure of 3C 345; the 
VERA array is sensitive to source structure because of its 
longest baselines in the JVN array. The improvements of 
the other parameters (e.g., peak intensity) arc the proof 
of efScacy of more-sophisticated calibration. 
3.3.3. Position shifts 

We evaluate here an astrometric capability from posi- 
tional reproducibility between the two epochs. Because of 
only a month, intrinsic changes of source structures were 
presumably negligible. The measured position offsets from 
the phase-tracking center on the target images arc listed in 
column (8) and (9) of Table 2. The improvements of posi- 
tional reproducibity by the BPR and more-sophisticated 
calibration are visualized in Fig. 4. The position differ- 
ences between the 1st epoch and 2nd epoch of the target 
with all the available antennas were steadily reduced in 
our calibration steps: 360+240 micro-arcsecond (/las) be- 
tween the normal phase-referenced images, 120+57 ^as 
between the BPR images, and finally, 62+50 /xas between 
the images by more-sophisticated calibration. Globally, 
the emission peak of 3C 345 appeared at ^^ 1 mas north- 
west from the phase-tracking center. These global posi- 
tion shifts include several origins, which will be discussed 
in Section 4.2. 

4. Discussions 

4-.1. Phase-referencing capabilities 

The BPR substantially achieves phase-referencing be- 
tween the C2-T pair without fast-switching between 
them. In our tests, 3C 345 and NRAO 512 were sepa- 
rated by 0.48 deg, 4.3-times smaller than the separation 



angle of actual fast-switching between 3C 345 and DA 426. 
A long-term phase drift from atmospheric/geometric er- 
rors, therefore, was expected to be 4.3-times smaller 
(Beasley & Conway 1995). However, resultant dynamic 
ranges of BPR images have been improved only 3.1-3.8 
times better than those of normal phase-referenced im- 
ages. This discrepancy was probably caused by sev- 
eral calibration errors (Section 2.5). The phase drifts 
seem to be successfully tracked with the sampling of 
C2 every 40 minutes. However, this sampling frequency 
might be risky in some level, as seen in Fig. 3 the 
phase of NRAO 512 and 3C 345 are systematically dif- 
ferent at about 9.2h UT. It is advisable to pair-switch 
every about 20 minutes, or to apply an another ob- 
serving schedule with more frequent C2 scans, such as 
-CI - C2 - CI - T - CI - C2 - CI - T - C1-. ... The 
observation case of indeed very weak C2, which is origi- 
nally supposed in the BPR, will be reported in a future pa- 
per. The BPR would give a better quality than that of the 
normal phase-referencing, even with a weak (~10 mJy) 
C2 with a signal-to-noise ratio of 5 in several scans using 
typical arrays at centimeter bands. 

We briefly discuss the over-correction of the solution ta- 
ble 2 in the BPR method, by comparison of 2nd epoch im- 
ages with all the antennas. In the normal phase-referenced 
images both of 3C 345 and NRAO 512 (not shown), very 
similar patterns of distorted contours are seen: positive 
and negative contours reside mainly at lower-right and 
upper-left side corners in the images, respectively [(i) in 
Fig. 2]. On the other hand, an opposite trend appears in 
the contours of the BPR image of 3C 345 [(j) in Fig. 2]. 
This is an effect of over-correction of solution table 2, 
which was designed for 2.57-degree separation of the Cl- 
C2 pair rather than 2.09-degree separation of the Cl-T 
pair. The distortion trend almost disappears in the im- 
age by more-sophisticated calibration using the optimized 
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Fig. 3. Plot of visibility phases phasc-refcrcnccd by primary calibrator DA 426 on Yamaguchi32ni-Iriki20m baseline. Open squares 
and filled circles represent phases of target 3C 345 and NRAO 512, respectively. Visibilities were averaged for 5 mimutes in this plot. 



table [(k) in Fig. 2]. 

4-.2. Astrometric capabilities 

The BPR and the more-sophisticated calibration sup- 
port relative astrometry, if observers use identical phase- 
tracking centers every epochs (Section 2). The consis- 
tency of the measured positions between our two epochs 
(Section 3.3.3) indicated that relative astrometry presum- 
ably worked well. The position differences among the 
epochs arc consistent with phase noises (Table 2). 

Dramatic revisions of the target position are seen par- 
ticularly in declination (Fig. 4). This indicates that the 
long-term phase drifts were responsible for a large part 
of the declination offsets. Geometric errors should be 
small in the VERA's four telescopes, Kashima34m, and 
Tsukuba32m, whose coordinates are strictly maintained 
by geodetic observations. The correlator model in the 
Mitaka FX correlator seems to have somewhat prediction 
errors in zenith excess-path lengths at troposphere and/or 
ionosphere. Note that a-priori calibration for such atmo- 
spheric errors is usually performed in the VERA project. 

There remain significant offsets ^ 1 mas west and ~ 
300 /ias north in both the epochs, even with the more- 
sophisticated calibration. Since such large global offsets 
could not be caused by any atmospheric/geometric errors 
in principle, an accumulation of intrinsic differences be- 
tween actual emission's positions and phase-tracking cen- 
ters of the target and two calibrators should be mainly 
responsible for the global offsets. This kind of offsets can- 
not be avoided even with the ICRF radio sources and their 
cataloged positions, which have uncertainties of ~ 0.3 mas 
(Ma et al. 1998; Fey et al. 2004), not only in the BPR and 
the more-sophisticated calibration but also in the normal 
phase-referencing. DA 426 (CI) has mas-scale jets east- 
southeast; 3C 345 (T) has mas-scale jets west-northwest. 
Hence, the global offsets pointing toward west-northwest 
may be reasonable. 

A position difference of 62±50 ^as between the two 
epochs in the images by more-sophisticated calibration 
might originate in the slightly misaligned allocation of the 
three sources. Strictly, the optimized solution table 2 re- 
acted to the sky position of C2'. This position was sepa- 
rated by 2.1 arcmin from T. The atmospheric and geomet- 
ric errors from this separation could cause a ~ 10-20 ^as 
difference between the two epochs. If the three sources 
perfectly aligned, the apparent positions of two epochs 
might coincide with each other within an error bar. 



4-3. Implications of the methods 

Although both the BPR and the more-sophisticated cal- 
ibration aim to retaliate against long-term phase drifts by 
reducing a separation angle, they are based on slightly 
different tactics. The BPR makes a weak source play a 
role of calibrator C2. Since self-calibration is done on 
C2, array's focus shifts to C2, closer to T than CI. The 
more-sophisticated calibration estimates the solution of 
self-calibration on a hypothetical calibrator C2' from the 
solution on C2. Array's focus finally shifts to C2', closer 
to T than C2. 

The key of the BPR is the presence of suitable C2 at 
a position less half the distance to T than that to CI. 
Without such a suitable C2, a BPR image will be worse 
than that of normal phase referencing. However, even a 
distant C2 may help the more-sophisticated calibration, if 
the sources align straightly. 

The key of the more-sophisticated calibration is the lin- 
earity of the phases Ai^atmo, A0goo, and A^pos. The 
linearity around sources have been already assumed in 
the method of ATMCA (Fomalont & Kogan 2005), which 
uses strong calibrators around a target to determine the 
phase gradient. The best performance can be obtained by 
both the more-sophisticated calibration and the ATMCA 
when sources align straightly. Even in case of misaligned 
source allocation, the calibration of ATMCA can be done 
by estimation of a two-dimensional phase gradient us- 
ing three or more calibrators. Because of the same as- 
sumption, the more-sophisticated calibration also can do 
it. Array's focus can be established exactly at T, when 

CiCa^ + CiC/ , i.e. riaC^ -t-rigCTO^, = C^f . Such a 
calibration table is derived from the two solution tables 
produced by self-calibrations on C2 and C3, in addition 
to that on CI. Not only bright calibrators but also weak 
ones are usable as secondary calibrators, when the BPR 
is applied before the more-sophisticated calibration or the 
ATMCA. Alternatively, with only one secondary calibra- 
tor, a good calibration can be achieved if assuming that 
the phase gradient is in elevation direction because the 
residual phase error is mostly tropospheric origin. This is 
supported by an option of the ATMCA. In such a process, 
observers should separate A^atmo from the other terms 
that would interact with the determination of the phase 
gradient. Fomalont & Kogan (2005) recommend observ- 
ing for more than 3 hours to find the true positions of 
calibrators, in order to make A0pos ~ 0. The term A(/)goo, 
however, would be harmful in arrays with poorly position- 
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Fig. 4. Improvements of reproducibility of target's (3C 345) 
position. The origin (0,0) of this plot is the phase-tracking 
center of the target. Labels A, B, and C show three different 
methods to measure positions. A: normal phase-referencing, 
B: BPR, C: more-sophisticated calibration. Three different 
symbols show different epochs and arrays. Squares repre- 
sent Ist-epoch measurements with VERA antennas. Triangles 
represent 2nd-epoch measurements with VERA antennas. 
Circles represent 2nd-epoch measurements with all antennas. 
The most calibrated images show a position difference of only 
62±50 /jas between two epochs (filled symbols). 

The BPR can make very weak sources available as caU- 
brators. In normal phase-referencing a signal-to-noise ra- 
tio of more than ~5 in half of a fast-switching period 
is required to obtain phase, delay, delay-rate solutions, 
while in the BPR a signal-to-noise ratio of ^^5 to ob- 
tain only a phase solution in about a quarter of a pair- 
swapping period is acceptable. This means that we can 
actually use calibrators more than a few times weaker 
than conventional calibrators. Consequently, the BPR 
brings lots of benefits, for example: (1) observers will be 
given more opportunities to make phase-referencing suc- 
cessfully at short centimeter and millimeter bands (~22- 
86 GHz), where targets normally have little chance of 
being accompanied by bright calibrators with a small 
separation angle. (2) Phase-referencing is available even 
to poorly sensitive or poorly position-determined anten- 
nas such as a space- VLBI. (3) Astrometry can be made 
between two weak sources, one of which has been self- 



calibrated. (4) Observers can easily obtain idealized cali- 
brator arrangements on the sky plane for the subsequent 
more-sophisticated calibration or the ATMCA, because 
of much larger surface number density in the faint-source 
sky. 

Practically, prior phase-referencing observations will be 
needed to find suitable calibrators around the target. We 
suggest that candidates are selected from flat-spectrum 
sources in the catalogs of the FIRST survey and so on, like 
an approach of 'a VLBA survey of flat-spectrum FIRST 
sources (Ulvestad ct al. 1999).' 

The BPR and the subsequent more-sophisticated cali- 
bration arc universal-designed. There is no requirement 
of upgrades for hardware or software. All observers have 
to do is scheduling of fast-switching with less frequent 
pair-swapping, and data reduction with bigradient cali- 
bration processes, which can be done only with the AIPS 
and editting solution tables. Even with antennas whose 
positions have not been well-determined, or even with a 
correlator without precise atmospheric models, very weak 
targets can be detected because of a much smaller sepa- 
ration from C2 or C2' than that from CI. A-priori cali- 
bration by geodetic- like observations (e.g., Mioduszewski 
& Kogan 2004) should bring independent improvements. 
Therefore, the combination of such prior calibrations and 
our method should achieve extremely high-quality phase- 
referencing. 

5. Summaries 

The bigradient phase referencing (BPR) allows us to 
utilize weak calibrators, one of which may be located at 
a position very closer to a target. The subsequent more- 
sophisticated calibration makes array's focus to shift to 
a hypothetical point much closer to a target. Thanks 
to much smaller separation angle, the phase-referencing 
quality of a target image will be dramatically improved. 
We describe the theory of the method and expected er- 
rors. For relative astrometry, observer must use identi- 
cal phase-tracking centers every epoch and the first self- 
calibrations on the calibrators using the structure mod- 
els at the phase-tracking centers. Our demonstrative ob- 
servation tests with strong sources (3C 345, NRAO 512, 
and DA 426) have shown its capabilities. Image dynamic 
ranges have been dramatically improved by a factor of 
about six, compared to normal phase-referencing. An as- 
trometric reproducibility is 62±50 //as between our two 
epochs in the most calibrated cases. The observation case 
of a weak calibrator, which is originally supposed in the 
BPR. will be reported in a future paper. 
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